NACA TN 3622 


CASE FILE 

caPT 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 3622 


PRELIMINARY STUDY OF SOME FACTORS WfflCH AFFECT THE 

STALL-FLUTTER CHARACTERISTICS OF THIN WINGS 

By A. Gerald Rainey 

Langley Aeronautical Laboratory 
Langley Field, Va. 



Washington 
March 1956 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 3622 

PRELIMINARY STUDY OF SOME FACTORS WHICH AFFECT THE 
STALL-FLUTTER CHARACTERISTICS OF THIN WINGS^ 

By A. Gerald Rainey 

SUMMARY 


The results of an exploratory^ analytical^ and experimental study 
of some of the factors which might he of importance in the stall flutter 
of thin wings are presented. The factors considered were Mach number^ 
Reynolds number^ density^ aspect ratio^ sweepback^ structural damping^ 
location of torsion nodal line^ and concentrated tip weights. 

The importance of aerodynamic torsional damping in the stall flutter 
of thin wings has been demonstrated through comparison of regions of 
negative torsional damping and regions of flutter. 

The results of a series of experiments on a thin wing tested at 
various lengths indicate that compressibility definitely alters the 
stall-flutter characteristics of wings of this type and that the com- 
pressibility effects appear to depend on the aspect ratio. 

A brief study of the inertia effects of concentrated tip weights 
indicates that such effects can be important} however^ the large number 
of parameters involved makes it difficult to generalize the results. 

An approximate analysis is presented for obtaining an estimate of the 
stall-flutter characteristics of particular wing -weight configurations. 

Some of the other factors considered were found to be more or less 
significant} however^ all the factors will require further study if 
their effects are to be more fully interpreted. 


INTRODUCTION 


In general^ as the angle of attack of a wing is increased to values 
near the angle of stalls the flutter velocity is decreased to values 
much lower than that observed at angles near zero lift. Since the 
potential-flow theories of classical flutter do not indicate an effect 


^Supersedes recently declassified NACA Research Memorandum L52D08 
by A. Gerald Rainey^ 1952. 
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of lift coefficient, the decrease in flutter velocity is presumed to be 
associated with nonpotential flow and has come to be known as stall 
flutter even though the wing may not be completely stalled in the usual 
sense of the word. Near zero lift the type of flutter which occvirs is 
usually a coupled, bending-torsion type, whereas flutter encountered at 
the higher angles of attack is usually predominantly a single-degree - 
of -freedom flutter occurring in the torsional mode. 

Much information exists concerning the stall flutter of moderately 
thick wings at low speeds (refs. 1 to 6) and the effects of several 
parameters on the stall flutter of thin propeller blades have been 
presented in reference 7> While stall flutter has not been considered 
a serious problem in the design of conventional wings (previous to 1950) 
experience with thin propeller blades indicates that stall flutter may 
become an important consideration in the design of thin wings. Since 
designers are considering the use of thin wings for high-speed airplanes 
it was considered desirable to reexamine the stall-flutter problem. 

The purpose of this paper is to present the resiilts of an explora- 
tory, analytical, and experimental study of some of the factors which 
might be of importance in the stall flutter of thin wings. The factors 
considered are Mach number, Reynolds number, density, aspect ratio, 
sweepback, structural damping, location of torsion nodal line, and 
concentrated tip weights. The results obtained by variation of any 
particular parameter are somewhat sketchy since the purpose of this 
investigation is to search for the significant parameters rather than 
to define completely the effect of the significant parameters. 


SYMBOLS 


a 


^n 


b 


c 


dE 

dt 


speed of sound, ft/sec 

initial amplitude of oscillation, radians 
amplitude of oscillation after n cycles, radians 
semichord perpendicular to leading edge, ft 
torsion spring constant per unit span, ft-lb/ft-radian 
chord, ft (2b) 

time derivative of energy, ft-lb/sec 
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Soi 

la 

Ka 

L 

M 

% 


n 


spanvise torsional deflection function 
first natural bending frequency^ cps 

first natural torsion frequency^ cps 

experimental flutter frequency^ cps 
damping coefficient 

structural damping coefficient for torsion mode 

moment of inertia per unit span about the elastic axis^ 
slug -ft ft 

spring constant of damping apparatus^ ft-lb/radian 
wing lengthy ft 

twisting moment per unit span^ ft-lb/ft 

complex aerodynamic moment per unit span per radian deflec- 
tion ft-lb/ft-radian 

aerodynamic moment per radian in phase with velocity^ 
ft-lb/radian 

coefficient of aerodynamic moment per radian in phase with 
velocity 

number of cycles 



effective inertia parameter where r^^ is the effective radius 
of gyration and k is the effective mass density ratio 


R Reynolds number 

S wing area^ ft 2 

V velocity, ft/sec 

X distance along span, ft 

a angular deflection, radians 
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a 

angular velocity, radians/ sec 


a 

angular acceleration, radians/sec^ 


Ui 

angle of attack at root, degrees 


% 

amplitude of torsional oscillation, 

, radians 

A 

angle of sweepback, degrees 


P 

density of test medium, slugs/ft^ 


o> 

frequency of oscillation, radians/sec 


first natural torsional frequency. 

radians/] 


= 2nf, 




METHODS AND APPARATUS 


The experimental investigation consisted of essentially two phases: 
(1) conventional flutter tests with the addition of the angle of attack 
as a primary variable and (2) aerodynamic damping measurements with the 
angle of attack again being a primary variable. 


Wind Tunnel 

The experiments were conducted in the Langley k. 5-foot flutter 
research tunnel which is of the closed-throat single-return type 
employing both air and Freon 12 as a testing medium at pressures from 
one atmosphere down to 1/8 atmosphere. During the course of this 
investigation the test section of the 4.5-foot-diameter flutter tunnel 
was altered to a rectangular section 2 feet by 4 feet to accommodate 
other types of investigations. The effects of this change in test 
section have not been evaluated but they are believed to be small. 


Models 

Five thin wings were used for the experiments and are illustrated 
in figure 1. Wing ninnber 1 was originally 24 inches longj however, in 
order to investigate the effects of aspect ratio, this wing was shortened 
in steps of three inches and the various configurations are designated 
wing number 1-24, 1-21, I-I8, and so forth. When equipped with a tip 
weight the wing designation is followed by the letter W, such as 1-21W. 
The model designations, properties, and test conditions are listed in 
table I. 
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Wing number 1 was of solid alumimom-alloy construction and was 
used in studying the effects of Mach number^ aspect ratio^ and con- 
centrated weights on stall flutter. This wing also was used to measure 
the aerodynamic damping moment for torsional oscillations . A series of 
four wings (designated 2^ 3^ and 5) varying in sweep angle and 
aspect ratio were employed to study the stall-flutter characteristics 
of swept wings. These wings were constructed of wood laminations on 
a stainless steel plate and the length -to -chord ratio was kept constant 
at 2.7» 

Each of the five wings was equipped with resistance-wire strain 
gages arranged to obtain a record of the flutter frequency and to 
indicate the relative magnitude of the strains in the bending and 
torsion degrees of freedom at flutter. 


Flutter Tests 

During the flutter tests the models were mounted as cantilever 
beams on a heavy^ rigid mount which could be rotated and clamped at 
any desired angle of attack. With the wing clamped at the desired 
angle of attack^ the airspeed of the tunnel was slowly increased -until 
a steady or divergent oscillation of the model occurred. At this point 
the tunnel conditions were observed and an oscillographic record of 
the strain-gage outputs was obtained. The velocity at this point was con- 
sidered to be the minimum flutter velocity for the wing tested at that 
particular angle of attack. In some cases the amplitudes encountered 
at this velocity were sufficiently mild so that the velocity could be 
increased further until some limiting condition was reached^ such as 
excessive vibration amplitudes^ excessive static bending load^ or 
maximum tunnel velocity. The same procedure was then repeated for 
angle-of-attack increments of 2^ or less at various angles throughout 
the range from 0*^ to . 

The effects of density^ Reynolds number^ and Mach nmber were 
studied by performing the above experiments at various pressures from 
atmospheric down to I/8 atmosphere in both air and Freon 12. 

The effects of aspect ratio and Mach number were investigated by 
shortening wing number 1 in steps of three inches. 

The effects of concentrated tip weights were studied, briefly, by 
clamping a solid ellipsoid, cast of an alloy of bismuth, tin, and lead, 

10 inches long and 0.75 inch in diameter to the tip of wings number 1-21 
and I-I8. The center of gravity of the weight was placed at the 
50-percent-chord line, 0.575 inch inboard from the tip in the plane of 
the wing. The characteristics of the wings equipped with the weight 
are listed in table I under the designations 1-21W and l-l8W. 
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Aerodynamic Damping Measurements 

The aerodynamic torsional damping moment was measured by the method 
of decaying oscillations (ref. 8). 

A special wing rnoimt (shown schematically in fig. 2) was designed 
and built for wing number 1-24 to allow torsional oscillations but 
restrict the bending mode. The root of the wing was restrained in 
torsion only by two relatively soft coil springs. The first natiiral 
torsional frequency of the wing-shaft -spring system was about 1/8 of the 
first natural torsional frequency of wing number 1-24 when mounted as a 
cantilever. In the special mount the wing tip was restrained in bending 
by a steel cable and, consequently, the first natural bending frequency 
of the wing in the damping mount was appreciably higher than in the 
cantilever condition. For these reasons the wing was considered as a 
rigid body (constant angle of attack spanwise) when oscillating at the 
first natural torsional frequency of the wing -shaft -spring system. 


The damping of the system was obtained by deflecting the wing- 
shaft-spring system in torsion to an amplitude of approximately 3° from 
the mean angle of attack. The wing was then released and the resulting 
oscillation was recorded. The damping coefficient was computed from 
the following relationship (ref. 9): 


1 , ^o 


and the aerodynamic damping moment was computed from 
Hi - (3-total “ Stare) Hx ~ Saero^ 

where gtotal damping coefficient obtained with air flow and 

Stare damping coefficient obtained with no air flow. The non- 

dimensional aerodynamic damping moment coefficient is defined as 


^ Py2sc 
2 

The aerodynamic torsional damping moment coefficient was obtained 
at various angles of attack by rotation of the entire wing -shaft -spring 
system to the desired angle of attack. The coefficient was determined 
as a function of v/bo3 for angles of attack of 0°, 4°, 8°, 12'^, l6°, 
and 20°. 
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The effects of the location of the axis of rotation or torsion 
nodal line were studied hy clamping the wing to the shaft so that the 
axis of the shaft coincided with the desired chordwise location. The 
nodal line positions studied were 25^ 32.5^ 50^ and 75 percent of the 
chord rearward of the leading edge. 

The maximum Mach number reached in the damping measurements was 
about 0 . 2 ; the maximum Reynolds number was about 0.8 X 10 ^. 


RESULTS Am DISCUSSION 

General Remarks Concerning the Effects of Lift 
Coefficient on the Flutter of Wings 


In general^ it has been found that as the angle of maximum lift is 
approached the flutter velocity is reduced appreciably from the value 
obtained at angles near zero lift. In addition^ the mode of flutter at 
the higher angles of attack is usually predOxminantly a single-degree- 
of-freedom flutter occirrring in the torsional mode. Since this type of 
flutter is associated with nonpotential flow^ the phenomenon^ as pre- 
viously noted^ has come to be known as stall flutter. 

A typical low Mach member variation of the flutter velocity coef- 
ficient with angle of attack for thin wings is shown in figure 3. The 
flutter velocity is plotted in the nondimens ional form v/ba^ as a 
function of the angle of attack at the root for wing number 2 . The 
lower boundaiy in figure 3 represents the velocity at which an oscilla- 
tion first started. The upper curve from = 6^ to = l8^ 
represents the velocity at which the oscillation reached dangerous 
amplitudes. The area labeled "flutter free” between 20 ^ and is 
bounded by the velocity coefficients at which the oscillation stopped. 

In this range of angle of attack the velocity was increased until the 
static bending load on the wing became excessive. The limiting velocity 
is indicated by the upper dashed curve. 

As the angle of attack is increased the flutter velocity coefficient 
decreases rapidly and reaches a minimum value near = l6^. The 
flutter velocity coefficient V/bc% at this minimum point is approxi- 
mately equal to 1 . 0 . The velocities listed in table I correspond to 
this minimum value for the various configurations and test conditions 
shown. Experimental evidence obtained on numerous propeller blades in 
reference 7 and on several thin wings in the present investigation 
indicates that the "rule-of-thumb^" v/bo^cc = 1.0^ for estimating the 
minimum flutter velocity of simple thin wings may be quite useful. Of 
course^ the exceptions to this elementary rule are numerous as will be 
indicated in the subsequent discussion of the various factors studied. 
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The remainder of the data shown in figure 3^ other than the minimum 
pointy do not lend themselves to generalities other than^ perhaps^ a 
few speculative remarks regarding the possible causes of the flutter. 

At zero angle of attack the flutter was a coupled bending-torsion type 
and the flutter frequency was at some value between the natural first 
bending frequency and the natural first torsional frequency. As the 
angle of attack was increased, the flutter-velocity coefficient decreased 
rapidly and the flutter frequency gradually approached the torsional 
frequency. This decrease in flutter velocity is probably due to a 
decrease in aerodynamic damping resulting from partial separation and 
boundary -layer time-lag effects. McCullough and Gault in reference 10 
indicate that, for steady flow, relatively sharp-edged airfoils exhibit 
a measurable region of separated flow near the leading edge at angles 
of attack as low as a = 1^. As the angle of attack is increased this 
separated-flow region is extended and, presumably, the boundai*y- layer 
time-lag effects become more pronoimced. The mechanism of this decrease 
in damping is not clearly understood and would be difficult to treat 
analytically. 

As the angle of attack is increased to values well beyond the 
steady-state stall angle (aj[ = 20^ to 24^), the vibration encountered 
resembles a forced vibration rather than a self-excited vibration. The 
amplitude of vibration increased gradually as the velocity was increased 
beyond the initial vibration point then, after reaching a maximum value, 
decreased to zero at the velocities indicated in figure 3- In this 
region of angle of attack, the vortex frequency, estimated on the basis 
of the projected wing area, was of the same order as the vibration 
frequency which was within a few percent of the natural torsional 
frequency. Although the motion under these conditions closely resembles 
a forced oscillation, it is difficult to separate the phenomenon at very 
high angles from the phenomenon at moderate angles and for that reason 
the oscillations at all angles of attack will be referred to as flutter. 


Some Considerations of Aerodynamic Damping 

As long ago as 1928 Frazer and Duncan (ref. l) hypothesized that 
the decrease in flutter speed with increasing angle of attack may be 
accounted for by a decrease of the aerodynamic torsional damping. In 
1943 Mary Victory (ref. 5) used a conventional flutter analysis altering 
only the torsional damping coefficients to calculate with good agree- 
ment several stall -flutter conditions obtained experimentally by Studer 
(ref. 2 ). The damping coefficients used in Victory*s analysis were 
obtained experimentally on a similar wing at about the same values of 
Reynolds number and frequency as the flutter data of Studer. 

Some of the aerodynamic torsional damping coefficients measured 
using wing number 1-24 are shown in figure 4(a). Contours of constant 
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values of the coefficient are shown in figure 4(a) as functions of 
v/bo) and angle of attack with the axis of rotation at the midchord. 

The curves above the zero damping line represent Increasing values of 
negative damping, except at the higher angles of attack where the damping 
decreased to zero then became positive again at higher values of V/bo). 
The area enclosed by the zero damping curves represents a region where 
a single-degree-of-freedom torsional instability is possible. 

For comparison, the actual flutter regions obtained for the same 
wing (no. 1-24) while mounted as a cantilever are shown in figure 4(b). 
One interesting feature of these data lies in the similarity of the 
flutter regions of figure 4(b) to the regions of negative damping shown 
in figure 4(a). This similarity is rather striking when it is realized 
that the damping data were obtained when the wing was mo'unted in a 
mechanism which provided appreciably different boundary conditions from 
those provided by the simple cantilever mount. The frequency of the 
wing mounted in the damping mechanism was about one-eighth of the 
torsional frequency of the wing moimted as a cantilever. The damping 
mechanism severely restricted the bending degree of freedom, whereas 
the cantilever mount, of course, allowed this degree of freedom. Tlie 
maximum Mach number of the damping meas-urements was about 0.2 as com- 
pared with a Mach number of about 0.75 for the flutter tests. The 
dissimilarity between the two sets of data at the lower angles of 
attack is attributed to compressibility effects which will be discussed 
later. In view of the similarity of the regions which occurred in spite 
of the widely different conditions which applied to the two determina- 
tions, it appears that the torsional aerodynamic damping is very impor- 
tant in the stall flutter of thin wings. 


Stall-Flutter Analysis 

The problem of calculating stall-flutter velocities has been a 
subject of much interest and some of the more promising approaches are 
discussed in reference 6. At the present time, no satisfactory stall- 
flutter analysis has been developed. This lack of an adequate analytical 
approach is due, primarily, to the difficulty in treating the air forces. 
The problem is somewhat similar to that of predicting the maximum lift 
coefficient. Presumably, if the oscillatory air forces were known for 
the proper range of significant variables, the analysis of stall flutter 
would be relatively simple. Unfortunately, the large number of signifi- 
cant variables makes it impractical to obtain adequate air-force coef- 
ficients except over very limited ranges. 

Observations of the stall-flutter characteristics of the unswept 
wings used in the present investigation indicate the possibility of 
applying a relatively simple analysis employing the measured aerodynamic 
torsional damping coefficients. The observations are: (l) The flutter 
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occurs in almost a pure "torsional mode wi"th very li"ttle transla'tion of 
the elastic axis and (2) the flutter frequency is very nearly equal to 
the first natural torsional frequency. (See table I.) Based on these 
observations a simplified single-degree-of-freedom analysis (developed 
in the appendix) yields the following relation: 



V \ 
n lbo)J 


= 0 


Since the structural damping is always positive, at flutter the aero- 
dynamic torsional damping must be negative and at least equal to the 
struct"ural torsional damping. 


If the aerodynamic torsional damping moment coefficients are 
available for the proper range of significant variables, the flutter 
velocity can be calculated. The aerodynamic damping moment coef- 
ficients measured in the present investigation are presented in fig- 
"ure 5 for the various angles of attack and axes of rotation. The 


2nu/ V / 

figure shows the val'ues of the quantity 1- — ) (the form most 

jt \^b(%/ 

venient for the calculations) as functions of v/bcD. For a given 
the quantity 


con- 


wing 



= Sa 



la [fa(x[] 


2 


dx 



2 

dx 


may be estimated, calculated, or measured, and the flutter velocity 
may be determined from the figures at each angle of attack where the 
aerodynamic damping is negative. Of course, at the lower angles of 
attack where the flutter can no longer be represented as a single- 
degree-of-freedom torsional vibration, the analysis is meaningless. 

Effects of structural damping .- One of the parameters which can 
cause deviations from the "rule of thumb" mentioned previously, that 
is, (v/batL)min = 1-0^ structural damping coefficient g^. This 

effect of structural damping (based on the preceding analysis) is 
illustrated in figure 6 where the minimum flutter velocity is shown 
as a function of the effective inertia parameter for several values 
of structural damping. The analytical results shown are for the axis 
of rotation at the midchord. The horizontal line for zero structural 
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damping indicates that the minimum flutter -velocity coefficient woTild 
he independent of the inertia parameter for this hypothetical condition. 
The addition of a small amo-unt of damping Increases the minimum flutter 
velocity coefficient appreciably and makes the coefficient a function 


of yir'y ■ ^ thin solid aluminum-alloy wing at sea-level density 


would have a value of 


near 20. Increasing the structural 


damping coefficient from 0.001 to 0.010 would approximately double the 
minimum flutter velocity coefficient. 


Effects of fluid density . - Variation of the density can cause 
deviations from the nominal minimum flutter velocity coefficient of 1.0. 
This effect is illustrated in figure 7, where the minimum flutter velo- 


city coefficient is plotted as a fvinction of the quantity 

which is Inversely proportional to density. The solid curve^ which 
was calculated for the midchord axis of rotation by the analysis 

previously discussed, rises sharply at values of representing 

high density then tends to level off at values corresponding to reduced 
density. The hypothetical thin aluminum-alloy wing mentioned earlier 



would have a value of 



of about 0.02 for a damping coefficient 


of g(x = 0.001 at sea level. The analysis indicates that the minimum 

flutter velocity coefficient would be increased from about 0.8 at sea 
level to 1.2 by decreasing the density to that at 40,000 feet altitude. 


Effects of location of the axis of rotation . - The effects of the 
location of the axis of rotation can be estimated with the aid of the 
aerodynamic torsional damping coefficients for various axes of rotation. 
Examination of figure 5 indicates a definite reduction in the extent of 
the negative damping regions for the forward locations of the axis of 
rotation. It appears, therefore, that the stall-flutter problem for 
wings having a torsion nodal line well forward may be greatly relieved 
as compared with the problem for wings having rearward locations of the 
torsion nodal line. 


Effects of Reynolds Number 

The effects of Reynolds number on stall flutter have not been 
thoroughly investigated} however, a few speculative remarks are possible. 
At low Mach numbers the aerodynamic characteristics of thin unswept wings 


12 


MCA TN 5622 


at high angles of attack are relatively insensitive to changes in 
Reynolds number. For this reason, it seems that Reynolds number might 
be relatively imimportant in the stall flutter of thin unswept wings 
at low Mach numbers. This supposition is at least partially verified 
by some of the data shown in figure T, which shows that the minimum 
flutter-speed coefficient for wing nimiber 1-24 increased from 


= 0.750 to 


= 0.846 over a range of density. The corre- 


sponding range of Reynolds number was from 3 -41 X 10^ to O .085 x 10^. 
The increase in flutter velocity coefficient with decreased density is 
in qualitative agreement with the trend predicted by the analysis. 


The effects of Reynolds number at high Mach ninnbers and for swept 
wings are appreciably more obscure and require firrther investigation. 


Effects of Mach Ntunber and Aspect Ratio on the Stall Flutter 

of Unswept Wings 

The "rule-of-thumb" mentioned earlier, that is, (r-^\ = 1.0 

\b<qx/ jjj 

(V)min = b(%, could be used to predict the minimum Mach nimiber at 
which flutter might occur if compressibility did not alter this minimum 
flutter velocity coefficient. The minimum flutter Mach number (v/a)jjj^^ 
would simply be bo^/aj however. Baker (ref. 7) found that, for thin 
propeller blades, (v/bu)a)niin function of Mach number. By varying 

the speed of sound a, it was shown that the minimum flutter Mach number 
became a nonlinear function of bc%,/a at the higher Mach numbers, and 
at sufficiently high values of bo^/a (approx. 0 . 5 ) no flutter occiarred 
up to the limit of the tests (M^^p = 1 . 4 ). 

The resiilts of the present investigation of compressibility effects 
are shown in figure 8 . The minimum flutter Mach number is shown as a 
function of the quantity bo^/a for wings number 1-24 to 1-15 . Since 
it was not possible to obtain a sufficiently wide variation of hcx^a 
by varying the speed of sound, it was necessary to increase 
by reducing the length of the original wing number 1-24. Consequently, 
the effects obtained may be thoiight of as combined effects of Mach 
rannber and aspect ratio. The results obtained on wing number 1-24 show 
a fairly linear relationship between the minimum flutter Mach niimber 
and bc^j/a and give no indication of an approaching limiting value 
of bcqx/a beyond which flutter would not occiur. The value of b<%/a 

reached in these tests was about 0.45, which is slightly less than the 
value of 0.5 presented in reference 7 as the value required to eliminate 
stall flutter of propeller blades. Wing munber 1-21 was tested at values 


NACA TN 5622 


13 


of bo^/a up to 0.5 and, although flutter was not eliminated at any 
value of boQx/a, the shape of the curve of minimum flutter Mach number 
as a function of bo^x/a (shown in fig. 8) indicated that a limiting 

value was approached at = 0.5* With wing number I-I 8 a new curve 

Si 


was defined with an apparent limiting value somewhere near = 0 . 59 * 

Si 

bow 

When wing number 1-15 was tested at = 0.7, no flutter was obtained 

a ’ 

up to the choking Mach number of the tunnel (M S 0.8); however, inter- 


mittent flutter did occur at — - = O. 67 , indicating that a more nearly 

a 

complete CTorve for this wing might have a shape similar to the dashed 
curve shown. 


The results of this series of experiments on one thin wing at 
various lengths indicate that compressibility definitely alters the 
stall-flutter characteristics of wings of this type. Furthermore, the 
compressibility effects appear to depend, to a large extent, on the 
aspect ratio. 


Effects of Sweepback 

The available data on the effects of sweepback on stall flutter of 
thin wings are limited to a series of experiments on wings number 2 , 3 ^ 

4, and 5- In. figure 9 the variation of the flutter velocity coefficient 
with angle of attack at atmospheric pressure is presented for the four 
wings varying in angle of sweepback from 0 ° to . 

The curves for the 0°, 15°^ and 30° wings are very similar, indi- 
cating no significant effect of sweep angle up to 30 ° for this series 
of wings. The 45° swept wing, however, did not experience flutter at high 
angles of attack even though the velocity was increased to about 80 per- 
cent of the flutter velocity obtained at zero angle of attack. The 
tests on the 45° wing were stopped at conditions indicated by the dashed 
curve because of the excessive static bending load on the wing. 

The results of the present experiments indicate the possible 
importance of sweepback; however, additional research will be required 
before the sweepback effects can be applied to othe'r configurations and 
test conditions. 
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Inertia Effects of Concentrated Tip Weights 


The inertia effects of concentrated tip weights on the stall- 
flutter characteristics of thin wings were investigated by equipping 
wings number 1-21 and I-I8 with a dense ellipsoid. High-density 
material was used to obtain a small volvime so that the aerodynamic 
effects of the mass would be a minimiam. The aerodynamic effects of 
large bodies on stall flutter have not been investigated. The weight 
reduced the natural torsional frequency of both wings approximately 
by a factor of 0.5* 

The variation of the flutter velocity coefficient with angle of 
attack for wings number 1-21 and I-I8 both with and without the weight 
is shown in figxire 10 . When tested without the weight both wings had 
a minimum flutter velocity coefficient of about 1.0, and this minimum 
occurred between = l6° and = 20 °. When tested with the weight, 

the minimum flutter velocity coefficient was increased to about 3-0 and 
the angle of attack at which this minimum occurred was decreased to 
about tti = 8° or 10 °. For comparison, the calculated variation (see 

section entitled "Stall-Flutter Analysis") for both wings with and 
without the weight are shown in the same figure. The calculated res\alts 
show the same trend as the experimental results in that the calculations 
indicate a large increase in the minimum flutter velocity coefficient 
and a decrease of the angle of attack at which this minimum occurs when 
the wing is equipped with a weight. 

The results shown in figure 10 indicate that the minimtun flutter 
velocity coefficient for these two wings was approximately tripled by 
equipping the wings with a concentrated tip weight, and the minimvim 
flutter velocity was approximately doubled. It should not be concluded, 
however, that the addition of a weight always increases the minimiun 
flutter velocity. The analysis on which the calculated results are 
based indicates that the minimum flutter velocity for a thin vinswept 
wing at low Mach numbers depends on the amoiint and distribution of mass 
or inertia, the fluid density, and the structural damping. For example, 
when wing number l-l8W was tested in Freon (high density) the minlmiom 
flutter velocity coefficient was about 1.5 as compared with 1.0 for wing 
number I-I8. The actual velocity, however, was a few percent lower for 
the weighted condition. 

Because of the many variables of importance in the stall flutter 
of wings with concentrated tip weights, it is not possible to generalize 
the results of this brief investigation. The Inertia effects of a 
particular wing-weight configuration can be estimated by use of figure 7 
where the minimum flutter velocity coefficient is shown as a function 


of the quantity 



The symbols represent experimentally 
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determined values for several wing -weight -density comhinations and the 
solid curve represents calculated minimum values for wings having a 
torsion nodal line near the midchord. The analysis agrees qualitatively 
with the experiments. 


CONCLUDING REMARKS 


The importance of aerodynamic torsional damping in the stall 
flutter of thin wings has been demonstrated through comparison of 
regions of negative torsional damping and regions of flutter. 

The results of a series of experiments on a thin wing tested at 
various lengths indicate that compressibility definitely alters the 
stall-flutter characteristics of wings of this type and that the com- 
pressibility effects appear to depend on the aspect ratio. 

A brief study of the inertia effects of concentrated tip weights 
indicates that such effects can be important; however^ the large number 
of parameters involved makes it difficult to generalize the results. 

An approximate analysis is presented for obtaining an estimate of the 
stall-flutter characteristics of particular wing -weight configurations. 

Some of the other factors considered (Reynolds number^ sweepback^ 
fluid density^ structural damping^ and location of torsion nodal line) 
were foimd to be more or less significant; however^ all the factors 
would require further study if their effects are to be more fully 
interpreted. 


Langley Aeronautical Labor at oiy 

National Advisory Committee for Aeronautics 
Langley Fields Va., March 3I, I952. 
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APPENDIX 


A simplified single-degree-of-freedom-flutter analysis using 
experimentally determined aerodynamic damping coefficients has been 
developed in order to provide some means for estimation of stall- 
flutter velocities of thin wings. 

The basic assumptions of the analysis are as follows: 

1. The motion of the wing is restricted to a single degree of 
freedom^ namely^ torsion. 

2. The flutter frequency o) is assumed to be equal to the natirral 
torsional frequency o^. 

If the approach used in reference 11 is followed^ the total moment 
contributed by an element of length dx is 


M = laa + Ca(l + iga)o- + ^aO. 


( 1 ) 


where Ma is the complex aerodynamic moment per unit span per radian 
deflection. 

The time rate of change of energy of each of the elements is equal 
to its total moment multiplied by its angular velocity. The rate of 
change of total energy for the entire wing having only torsional motion 
is simply the sum of the energy rates contributed by each of the 
elements 


dt 


"L 

Md. dx 


0 


( 2 ) 


At the flutter condition the time rate of change of energy must equal 
zero. Substituting the expression for the moment, equation (l), into 
the expression for energy, equation (2), and setting the total rate of 
change of energy equal to zero produces 



[la.a + Ca(l + iga)°- + Ma<0d dx = 0 


dt 


( 3 ) 
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The motion of each element is harmonic and can be expressed as 


a 


= 




where represents the spanwise mode shape for torsional vibra- 

tions. It is assumed that the aerodynamic terms do not appreciably 
change the mode shape. 

If this relation is used and Ccc = la^^ is defined^ equation (3) 
becomes 



0^2(1 + ig^)J 



la |£a(x)] ^dx + 





Equation (4) can be separated into a real and an imaginary equation 
both of which must vanish to provide a solution to (^). Since it has 
been observed that in most cases of stall flutter o) = the real 
part of equation (4) must approach zero. Normally the real part of 
equation (4) is used to provide the flutter frequency but^ since the 
flutter frequency has been assumed to be equal to the natural torsional 
frequency^ the solution of the real part of equation (4) will not be 
required. 


The imaginary part of equation (4)^ solutions of which provide the 
flutter velocity^ is 




PL 


la [pa( xl] 


"dx + 




(imag. )Ma [fa(x]] ^dx = 0 (5) 


or 


Ia[fa(x)]^dx + 2py2 |fa(x3 ^dx = 0 (6) 

Jo Jo 


The first part of equation (6) represents the structural damping 
energy whereas the second part represents the aerodynamic damping energy. 
If it is assvuned that a representative value of the semichord can be 
chosen for purposes of defining that is, not a function 
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of 


or 


where 


then equation (6) can be simplified as 



represents the effective or mean value of the important 


k: 



inertia parameter 
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TABLE I 
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MODEL PARAMETERS AND CONDITIONS AT TBE MINIMUM FLUTTER VELOCITY 


Model 

designa- 

tion 

Aspect 

ratio 

A 

(deg) 

fhl 

(cps) 

fai 

(cps) 

Sa 

b 

(ft) 

“•i 

(deg) 

V 

( ft/ sec) 

v/Ht 

Mach 

number 

R X 10"^ 

a 

(ft/sec) 

bo^/a 

Percent 

Freon 

Density 

P 

( slue/ 

cu ft) 

f 

•^exp 

(cps) 

Test 

section 

(ft) 

1-24 

6.0 

0 

17.3 

102.3 

0.0010 

0.333 

20 

171.8 

0.803 

0.149 

0.651 

1152.8 

0.186 

0 

0.00222 

100.2 

4.5 diam. 

1-24 

6.0 

0 

17.3 

102.3 

.0010 

.333 

22 

173.5 

.811 

.151 

.332 

1149.1 

.186 

0 

.00112 

102.2 

4.5 diam. 

1-24 

6.0 

0 

17.3 

102.3 

.0010 

.333 

20 

178.1 

.832 

.154 

.162 

1157.9 

.185 

0 

.00054 

103.0 

4.5 diam. 

1-24 

6.0 

0 

17.3 

102.3 

.0010 

.333 

20 

181.1 

.846 

.160 

.085 

1135.5 

.189 

0 

.00027 

103.0 

4.5 diam. 

1-24 

6.0 

0 

17.3 

102.3 

.0010 

.333 

20 

160.6 

.750 

.322 

3.411 

499.0 

.429 

97.0 

.00846 

92.1 

4.5 diam. 

1-24 

6.0 

0 

17.3 

102.3 

.0010 

.333 

20 

170.5 

.796 

.336 

1.853 

507.8 

.422 

94.8 

.00441 

97.4 

4.5 diam. 

1-24 

6.0 

0 

17.3 

102.3 

.0010 

.333 

20 

182.6 

.853 

.344 

.849 

531.0 

.403 

87.1 

.00198 

100.9 

4.5 diam. 

1-24 

6.0 

0 

17.3 

102.3 

.0010 

.333 

20 

188.4 

.880 

.326 

.346 

578.0 

.370 

69.3 

.00083 


4.5 diam. 

1-21 

5.25 

0 

22.0 

121.4 

.0010 

.333 

20 

217.3 

.855 

.215 

.137 

1010.0 

.252 

6.7 

.00035 

120.0 

2x4 

1-21 

5.25 

0 

22.0 

121.4 

.0010 

.333 

19 

272.3 

1.071 

.537 

5.341 

507.0 

.502 

93.9 

.00781 

115.8 

2x4 

1-21 

5.25 

0 

22.0 

121.4 

.0010 

•333 

18.5 

291.1 

1.145 

.568 

2.846 

513.0 

.496 

90.7 

.00395 

116.7 

2x4 

1-21 

5.25 

0 

22.0 

121.4 

.0010 

.333 

20 

217.3 

.855 

.392 

.911 

555.0 

.458 

81.1 

.00184 

120.7 

2x4 

1-21W 

5.25 

0 

— 

68.0 

.0057 

.333 

9.5 

358.1 

2.513 

.324 

.204 

1107.0 

.128 

2.2 

.00032 

67.7 

2x4 

1-18 

4.5 

0 

29.2 

143.7 

.0012 

.333 

18 

287.9 

.956 

.251 

.120 

1147.0 

.262 

0 

.00024 

143.7 

2x4 

l-l8 

4.5 

0 

29.2 

143.7 

.0012 

.333 

18.3 

299.6 

.996 

.600 

2.926 

514.5 

.585 

92.3 

.00385 

l4o.o 

2x4 

1-18 

4.5 

0 

29.2 

143.7 

.0012 

.333 

18.3 

290.0 

.964 

.561 

2.611 

517.0 

.582 

91.6 

.00365 

l4o.o 

2x4 

1-18 

4.5 

0 

29.2 

143.7 

.0012 

.333 

18.3 

287.3 

.955 

.538 

1.114 

534.0 

.564 

85.2 

.00163 

145.0 

2x4 

1-18 

4.5 

0 

29.2 

143.7 

.0012 

.333 

20.2 

268.9 

.894 

.481 

.916 

559.0 

.538 

7^.5 

.00149 

145.0 

2x4 

l-l8w 

4.5 

0 

17.4 

81.9 

.0052 

.333 

8.0 

568.5 

3.275 

.500 

.269 

1137.0 

.151 

0 

.00027 

81.3 

2x4 

l-l8w 

4.5 

0 

17.4 

81.9 

.0052 

.333 

14.0 

262.3 

1.510 

.476 

.102 

551.0 

.311 

78.5 

.00172 

81.3 

2x4 

1-15 

3.75 

0 

4l.i 

179.1 

.001 

.333 

18.3 

391.7 

1.043 

.697 

1.420 

562.0 

.668 

77.9 

.00147 

175.0 

2x4 

1-15 

3.75 

0 

4l.i 

179.1 

.001 

.333 

18.3 

427.2 

1.139 

.794 

1.308 

538.0 

.698 

85.2 

.00142 

175.0 

2x4 

2 

5.4 

0 

10.1 

36.2 

.017 

.208 

16.0 

48.0 

.973 

.042 

.116 

1146.0 

.04l 

0 

.00226 

35.7 

4.5 diam. 

2 

5.4 

0 

10.1 

36.2 

.017 

.208 

i4.o 

56.0 

1.184 

.049 

.068 

1147.0 

.04l 

0 

.00113 

35.2 

4.5 diam. 

2 

5.4 

0 

10.1 

36.2 

.017 

.208 

8.0 

109.9 

2.322 

.097 

.073 

1131.0 

.042 

0 

.00060 

35.4 

4.5 diam. 

3 

5.04 

15 

10.0 

40.3 

.02 

.208 

18.0 

51.1 

.967 

.045 

.134 

1134.0 

.047 

0 

.00231 

37.8 

4.5 diam. 

3 

5.04 

15 

10.0 

40.3 

.02 

.208 

l4.o 

54.1 

1.025 

.048 

.070 

1125.0 

.047 

0 

.00114 

37.9 

4.5 diam. 

3 

5.04 

15 

10.0 

40.3 

.02 

.208 

l4.o 

73.2 

1.388 

.065 

.049 

1129.0 

.047 

0 

.00058 

38.9 

4.5 diam. 

4 

4.05 

30 

10.6 

38.1 

.02 

.208 

20.0 

58.7 

1.177 

.052 

.173 

1132.0 

.044 

0 

.00233 


4.5 diam. 

4 

4.05 

30 

10.6 

38.1 

.02 

.208 

16.0 

70.9 

1.423 

.062 

.099 

1139.0 

.044 

0 

.00117 


4.5 diam. 

4 

4.05 

30 

10.6 

38.1 

.02 

.208 

8.0 

203.2 

4.079 

.179 

.152 

1136.0 

.044 

0 

.00059 

38.8 

4.5 diam. 

5 

2.7 

45 

12.4 

40.9 

.02 

.208 

8.0 

354.3 

6.616 

.311 

1.160 

ll4l 

.047 

0 

.00217 


4.5 diam. 
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Figure 1 .- Flutter models. 
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Coil springs 



Bending restraint 


Spring support frame 


Figure 2 .- Aerodynamic damping mechanism. 
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Angle of attack, 


Flgirre 3.- Flutter regions for wing number 2 at one atmosphere in air. 
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(a.) Con'touxs of const/ant aerodynamic 'torsional damping coefficien't as 
functions of V/lxs and aj^. Axis of rotation at 0.50c, 


V 

bu<. 



(b) Flutter regions for -wing number 1-24 
at 1/8 atmosphere in air. 

Figure 4,- Comparison of flutter regions and damping contours. 
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(a) Axis of rotation at 0.25c. 


Figtire 5»- Variation of 


Jt \ba)/ 


with — for various angles of attack. 
hCD 


2m. 
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(b) Axis of rotation at 0.325c. 
Figiure 5.- Continued. 
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(c) Axis of rotation at 0.50c. 


Figiire 5-- Continued 


2m . 
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ro 

cx» 
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bcj 


(d) Axis of rotation at 0.75c. 


Figure 5*- Concluded 
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Figure 6 .- Calculated variation of ininiim:iin flutter velocity coefficient 
with effective inertia parameter for various values of structural 
damping. Axis of rotation at 0 . 5 c. 


5 


o 


min 


Wing 

No. 

C 1-24 
□ 1-24 
A 1-21 
O 1-21W 
fc. 1-lS 

G 1-lSW 
A 1-l^W 


^ Calculated 



Figure 7.- Comparison of experimental and calculated minimum flutter 
velocity coefficients for wings number 1-24 to l-l8 at various 
densities with and without a concentrated tip weight. 
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XJ 

o 

C0 
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s 
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a 
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Wing Test Aspect 
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b<4. 


boDa 

Figure 8.- Variation of minimum flutter Mach niimher with . 

3. 
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V 





Angle of attack, OC 


Figure 9.- Variation of flutter velocity coefficient with angle of attack 
for various angles of sweepback at 1 atmosphere in air. 
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V 

boi« 



Angle of attack, 


Figure 10 .- Comparison of experimental and calculated flutter velocity- 
coefficients for wings number 1-21 and I-I8 with and without concen- 
trated tip weight at I/8 atmosphere in air. 
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